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Abstract 

Hopf solitons in the Skyrme-Faddeev model — S' 2 -valued fields on M 3 with Skyrme 
dynamics — are string-like topological solitons. In this Letter, we investigate the 
analogous lattice objects, for S' 2 -valued fields on the cubic lattice Z 3 with a nearest- 
neighbour interaction. For suitable choices of the interaction, topological solitons 
exist on the lattice. Their appearance is remarkably similar to that of their continuum 
counterparts, and they exhibit the same power-law relation E ~ ciL 3 / 4 between the 
energy E and the Hopf number H. 
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One of the simplest 3-dimensional systems admitting topological solitons is the 
0(3) sigma model — in other words, where the field is a unit 3-vector ip{x J ) on 
M 3 . Given the usual boundary condition ip{x J ) —> T (constant) as |x J | —> oo, such 
configurations are classified by their Hopf number H[ip] e -rr^^S 2 ) = Z. There is an 
integral formula for H[ip]: first define Fj k = ip ■ {djip ) x {d k ip)\ then find Aj such that 
Fjk = djAk — 8 k Aj (this is always possible); and finally, compute 



An energy functional of the form E = f[\dip \ 2 + V{ip)\d 3 x does not admit sta¬ 
tionary soliton solutions with nonzero H[ip ]: for example, E has no smooth critical 
points, and there are no stationary solutions of the Landau-Lifshitz equation [T]. 
But with a modified E, one can have interesting soliton solutions. The best-known 
example comes from adding a Skyrrne term nig, so that 

E[ip\ = J [{djip) ■ {djip) + i Fj k Fj k } d 3 x. (2) 

It is then believed that E has a smooth minimum in each topological class, with 
E ~ cH 3 / A for some constant c; and there is considerable evidence supporting this 
conjecture, both analytic HEJEI and numerical 0@EHE|. These solitons can be 
visualized as closed curves, which may be linked or knotted. In particular, if p is 
the point on S 2 corresponding to the boundary value T, and q G S 2 is the antipodal 
point, then the soliton may be viewed as the closed curve ip~ l {q) in M 3 , which links 
| H\ times around the open curve ip^ 1 {p). 

In this Letter, we consider a lattice rather than a continuum system: namely, 
the field consists of a unit vector cp defined at each point of the three-dimensional 
cubic lattice Z 3 . At first sight, one then loses all the topological features; but by 
restricting to “well-behaved” configurations, many of the topological properties can 
be maintained. There are many well-known examples of this, in other systems ad¬ 
mitting topological solitons. One such is the 0(3) model on the 2-dimensional lattice 
Z 2 : for a suitable set of configurations, one can define the topological charge un, 
and (depending on the choice of dynamics) one can have a Bogomolny bound on the 
energy, and the existence of stable topological solitons 02 Eg. We shall see below 
that analogous results hold for the 3-dimensional case. In particular, it makes sense 
to talk about Hopf solitons on Z 3 ; it turns out that they resemble their continuum 
counterparts, and that their energy has the same behaviour E ps ciL 3 ^ 4 as in the 
continuum case. 

So let us consider a spin (unit 3-vector) cp = 4>i,j,k defined at each site of the three- 
dimensional cubic lattice Z 3 . Equivalently, one may think of a map cp : Z 3 —► S 2 . 
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The boundary condition at infinity is (pi.j t k —>'£:= (0,0,1) as i 2 + j 2 + k 2 —> oo. We 
assume there to be an isotropic nearest-neighbour interaction, so the energy E[(p\ of 
a configuration has the form 


E i4>\ = \f(4>i+i,j,k 

i,j,k 



(3) 


where / is a suitable function. We may take /(1) = 0, so that the constant field 
(f>i j k = has zero energy; and /'(1) = —1, which sets the energy scale. 

The simplest such function is /(£) = 1 — which corresponds to the usual 
Heisenberg model; but for this choice of /, the only minimum of E[(f>\ is the constant 
field 4>(x 3 ) = <h. In order for interesting local minima of the energy to exist, we need 
to make it energetically more unfavourable for (p to become “discontinuous”, in the 
sense that neighbouring spins point in wildly different directions. For example, this 
can be achieved by taking /(£) = (1 — £) + a( 1 — £) 2 , with the parameter a being 
sufficiently large. 

Another way to view the situation is that we want to restrict to “continuous” 
configurations (p, for which a Hopf number H[(j)\ € Z can be unambiguously defined. 
This is a familiar idea: as mentioned above, for example, one can define the winding 
number of a spin-field on Z 2 by excluding certain exceptional configurations ED- In 
the present case, one may think of Z 3 as embedded in M 3 ; if cp : Z 3 —> S 2 can be 
extended to a continuous map ip : M 3 —> S 2 , in a way which is unambiguous up to 
homotopy, then H[cp] can be defined as H[ip\. In order for such an unambiguous 
interpolation ip to exist, we need a “continuity” condition on cp. Let us impose the 
following condition: that the angle between any pair of nearest-neighbour spins is 
acute. Continuous fields 0 are those which satisfy this condition: </>j+ i,j,k ' fiij,k > 0 
for all i,j, k, with a similar inequality for the j-links and the A:-links. 

To see what might happen if cf> becomes discontinuous, consider a face of the 
lattice, say P = {(0,0,0), (1, 0,0), (0,1,0), (1,1,0)}. Its image on S 2 is a spherical 
quadrilateral. If neighbouring spins were allowed to be orthogonal, then the four 
vertices {0o,o,o> 01,0,0) 0 o,i,Oj ^l.i.o} of this quadrilateral could lie on the equator of 
S 2 . But then it would be unclear how to interpolate: should the image of the interior 
of the face P be the northern hemisphere of S 2 or the southern hemisphere? For 
continuous fields this ambiguity cannot occur, and the interpolation ip is well-defined 
up to homotopy. In the two-dimensional case (maps from Z 2 to S 2 ), one can then 
define the winding number directly, by adding up the signed areas P(P) of the 
spherical quadrilaterals corresponding to each face mi- For the three-dimensional 
case, computation of H[cp\ is not quite so straightforward: we shall return to this 
below. The basic fact is that, just as in the continuum case, the space of continuous 


3 



configurations is disconnected, and its components are labelled by the Hopf number 
H[</>] G Z. 

As mentioned above, the simplest choice /(£) = 1 — £ of inter-spin potential does 
not permit the existence of static solutions (local minima of E) with non-zero Hopf 
number H . If one starts with a continuous configuration cj) having H\c/)\ A 0, and 
allows it to flow down the energy gradient, then cj) becomes discontinuous, and the 
topology is lost. If this is to be avoided, the function /(£) needs to contain higher 
powers of 1 — The simplest choice is 

/(£) = ( l-0 + «(l-0 2 - (4) 

The parameter a has to be large enough to avoid the instability referred to above, and 
numerical experiments indicate that a value of a = 25 is sufficient for this (whereas 
a = 16, say, is not). In what follows, we adopt (J3J) with a = 25. 

The system defined by @0 was investigated numerically: the procedure con¬ 
sisted of mimimizing the energy E , using a conjugate-gradient method. The quantity 
min{0 • cj) + } (the minimum being taken over all links) was monitored, to check that 
it remained positive — in other words, that <f> remained continuous as it flowed down 
the energy gradient. The computation was done for finite lattices of size N 3 , with 4> 
set equal to $ = (0, 0,1) at the boundary of the lattice, for a range of values of N 
up to N = 22; the results were then extrapolated to the unbounded ( N —> oo) limit. 

In addition, the Hopf number H\cj)\ was monitored. One can get an approximate 
value for H[4>\ by using a discrete version of the formula (0, as follows. The first step 
is to compute, for each face P of the lattice, the signed area E(P) of the corresponding 
spherical quadrilateral. This object J~, which assigns a number E{P) to each face P, 
is a 2-cocycle on Z 3 , and is the analogue of the 2-form F in 0 . Given E, one can 
compute a 1-cochain A satisfying 5A = E\ this object A assigns a number A(L) to 
each link L, and is the analogue of the 1-form A in 0. Then 

^° M = 12^ ( 5 ) 

P±L 

gives an approximation to H\(f)\. By interpolating (j> to a lattice with half the lattice 
spacing (ie. with eight times as many sites) and repeating the calculation, one gets a 
better approximation Hi[(j)]; and this procedure can be iterated to obtain a sequence 
H n which converges to H. In fact, the quantity H' = (4 Hi[(j>] — Hq[4>])/3 is already 
within 1% of the true value H, for the fields considered here; and it was H' which 
was monitored to check that the soliton had the required Hopf number. 

The initial configurations were taken to be axially-symmetric, with the symmetry 
axis chosen so as not to be aligned with the lattice. Axisymmetric fields [7j are la¬ 
belled by two integers (m, n), with mn = E[; for each value of H, various possibilities 
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for (m, n) were tried, and the one leading to the lowest-energy soliton selected. For 
example, in the H = 4 case, the lowest-energy configuration is of type (2,2), as in 
the continuum case m , and it retains the symmetry of the initial configuration. For 
H = 6, by contrast, the symmetry of the initial axisymmetric (3, 2) configuration is 
lost. In some cases (such as H = 2), the initially non-aligned soliton rotates during 
the minimization so as to become aligned with the lattice. In other cases (such as 
H = 1), the non-aligned soliton has a lower energy than the corresponding aligned 
one (which one can obtain by starting with an aligned initial configuration). The 
picture seems to be that there are even more local minima, in each topological class, 
than in the continuum case — with these additional solutions arising because of the 
anisotropy of the lattice. 

The results of the numerical experiments may be summarized as follows. For each 
value of the Hopf number 1 < H < 6, there exists at least one stable local minimum 
of E[<p\, and the minimum energies are plotted in Figure 1, on a log-log scale. The 
dashed line is the line through the H = 1 data point having slope 3/4; so it is clear 
that the same power law E ~ ciF' 4 / 4 holds as in the continuum Skyrme-Faddeev 
system. The corresponding configurations <j) are illustrated in Figure 2. Recall that 
the boundary condition is (f> —> (0,0,1) as r -> oo, and in the continuum we can 
visualize the soliton as being located on the closed curve </>3 = — 1 corresponding to 
the antipodal point on the target space S 2 . On the lattice, we obtain an analogous 
picture by plotting the surfaces = k for some suitable constant k. These surfaces 
are plotted in Figure 2, using k = —0.6. In these pictures, the lattice spacing is unity; 
so one sees that the solitons are typically spread over a width of 10-20 lattice units 
(depending on their shape and on the value of FI). It is remarkable that the pictures 
closely resemble the corresponding ones in the continuum case HEED!, despite the 
systems being quite different. 

These results are just for one particular choice of the inter-site potential, namely 
But various other choices of / have been investigated as well, and they lead to 
similar results: for example, the choice m 

OO 1 

/{{) = -log? = £-(l-{)\ (6) 

1 ^ 
n= 1 

which has the feature that E[(j)\ —» oo as </> approaches any discontinuous configu¬ 
ration. The conclusion, therefore, seems to be that solitons in 3-dimensional 0(3) 
models have “universal” features: provided they exist at all, their appearance, as well 
as the 3/4 power law for their energy, are the same for a wide variety of choices of 
dynamics. It might be interesting to investigate whether there are analytic reasons 
for generic features such as these. 
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Figure 1: The energy E of minimal-energy configurations, as a function of the Hopf number 
H, for 1 < H < 6. The dashed line has slope 3/4. 
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